Challenging Magnetic Field Dependence of the Residual Resistivity of the 

Heavy-Fermion Metal CeCoIn 5 
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An explanation of paradoxical behavior of the residual resistivity po of the heavy-fermion metal 
CeCoIns in magnetic fields and under pressure is developed. The source of this behavior is identified 
as a flattening of the single-particle spectrum, which exerts profound effects on the specific heat, 
thermal expansion coefficient, and magnetic susceptibility in the normal state, the specific heat 
jump at the point of superconducting phase transition, and other properties of strongly correlated 
electron systems in solids, ft is shown that application of a magnetic field or pressure to a system 
possessing a flat band leads to a strong suppression of po. Analysis of its measured thermodynamic 
and transport properties yields direct evidence for the presence of a flat band in CeCoIns. 
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Measurements of the resistivity p(T) in external 
magnetic fields H have revealed a diversity of low- 
temperature behaviors of this basic property in heavy- 
fermion (HF) metals, ranging from the familiar Landau 
Fermi liquid (LFL) character to challenging non- Fermi 
liquid (NFL) behavior— ~— . The resistivity p(T) is fre- 
quently approximated by the formula 



p(T) = Po + AT n , 



(1) 



where po is the residual resistivity and A is a T- 
independent coefficient. The index n takes the values 
2 and 1, respectively, for FL and NFL behaviors and 
1 < n < 2 in the crossover between. The term po is ordi- 
narily attributed to impurity scattering. The application 
of the weak magnetic field is known to produce a posi- 
tive classical contribution oc H 2 to p arising from orbital 
motion of carriers induced by the Lorentz force. How- 
ever, when considering spin-orbit coupling in disordered 
electron systems where electron motion is diffusive, the 
magnctoresistivity may have both positive (weak local- 
ization) and negative (weak anti-localization) signs^. 

Our focus here is on the compound CeCoIns, whose 
H — T phase diagram is drawn schematically in Fig. [T] 
Its resistivity p(T, H) differs from zero in the region be- 
yond the solid curve that separates the superconducting 
(SC) and normal states. Above the critical temperature 
T c of the SC phase transition, the zero-field resistivity 
p(T, H = 0) varies linearly with T. On the other hand, 
at T — >• and magnetic fields H > H C 2 ~ 5 T, the 
curve p(T,H C 2) is parabolic in shapoA^.. As studied ex- 
perimentally, CeCoIns is one of the purest heavy-fermion 
metals. Hence the applicable regime of electron motion 
is ballistic rather than diffusive, and both weak and anti- 
weak localization scenarios are irrelevant. Accordingly, 
one expects the -H-dependent correction to po to be pos- 
itive and small. Yet this is far from the case: specifi- 



cally p (H = 0) ~ 1.5 pflcm, while p (H = 6T) ~ 0.3 




Magnetic field H, arb. units 



FIG. 1: (color online). Schematic T — H phase diagram of 
CeCoIns. The vertical and horizontal arrows crossing the 
transition region marked by the thick jagged lines depict the 
LFL-NFL and NFL-LFL transitions at fixed H and T, respec- 
tively. As shown by the solid curve, at H < H C 2 the system is 
in its superconducting (SC) state, with H c o denoting a quan- 
tum critical point hidden beneath the SC dome where the 
flat band could exist at H < H c q. The hatched area with the 
solid curve T CIOSS (H) represents the crossover separating the 
domain of NFL behavior from the LFL domain. A part of 
the crossover marked with the dots is hidden in the SC state. 
The NFL state is characterized by the entropy excess S, of 
Eq. (El). 



To resolve this paradox, we suggest that the electron 
system of CeCoIns contains a flat band. Flattening 
of the single-particle spectrum e(p) is directly relevant 
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to the problem addressed since, due to Umklapp pro- 
cesses, quasiparticles of the flat band produce a contri- 
bution to po indistinguishable from that due to impu- 
rity scattering^. Furthermore, it is crucial that the flat 
band somehow becomes depleted at T — > and H = 6 
T, to avoid contradiction of the Nernst theorem. This 
depletion entails a dramatic suppression of the flat-band 
contribution to po- 

Before proceeding to the analysis of this suppression 
and its relevance to the behavior observed in CeCoIns, 
we call attention to salient aspects and consequences of 
the flattening of e(p) in strongly correlated Fermi sys- 
tems. The theoretical possibility of this phenomenon, 
also called swelling of the Fermi surface or fermion con- 
densation, was discovered two decades ago&~— (for recent 
reviews, seeii~— ). It has received new life in the concep- 
tual framework of topological matter, where it is charac- 
terized by nontrivial topology of the Green function in 
momentum space and associated with topologically pro- 
tected flat bands^~— . At T = 0, the ground state of a 
system having a fiat band is degenerate, and therefore 
the occupation numbers rt*(p) of single-particle states 
belonging to the flat band, which form a so-called fermion 
condensate (FC), are continuous functions of momentum 
that interpolate between the standard LFL values {0, 1}. 
This leads to an entropy excess 

5* = -^n»(p)lnri*(p) + (f-n*(p))ln(l-n*(p)), (2) 
p 

which does not contribute to the specific heat C(T). 
However, in contrast to the corresponding LFL en- 
tropy, which vanishes linearly as T — > 0, S*» pro- 
duces a T-indcpcndent thermal expansion coefficient a oc 
~dS*/dP—, where P is the pressure. In its normal state, 
CeCoIns does in fact exhibit a greatly enhanced and al- 
most T-independent thermal expansion coefficient 2 ^, so it 
is reasonable to assert that it possesses a flat band. Anal- 
ysis of the experimental data on magnetic oscillation o 21 ! 22 
supports this assertion. CeCoIns is found to have two 
main Fermi surfaces. The a-sheet is observable at all 
fields down to H C 2, while the magnetic oscillations asso- 
ciated with the /3-sheet become detectable only at mag- 
netic fields H > 15 T, behavior in accord with the posited 
flat character of this band. 

In the theory of fermion condensation, the aforemen- 
tioned ground-state degeneracy is lifted at any finite tem- 
perature, where FC acquires a small dispersion propor- 
tional to T, the spectrum being given by-i^ 

e (p,n,)=rin 1 ' 7 ;* (p) . (3) 

However, the lifting of the degeneracy does not change 
the FC occupation numbers n* (p) , implying that the en- 
tropy excess 5** would persist down to zero temperature. 
To avert a consequent violation of the Nernst theorem, 
FC must be completely eliminated at T — > 0. In the most 
natural scenario, this happens by means of a SC phase 



transition, in which FC is destroyed with the emergence 
of a pairing gap A in the single-particle spectrurn^^r— . 
We propose that this scenario is played out in CeCoIns at 
rather weak magnetic fields H <C H C 2 , providing for elim- 
ination of the flat portion in the spectrum e(p) and the re- 
moval of the entropy excess S 1 * . In stronger external mag- 
netic fields sufficient to terminate superconductivity in 
CcCoIn5, this route becomes ineffective, giving way to an 
alternative scenario involving a crossover from FC state 
to a state having a multi-connected Fermi surfac o 11 ^ 7 " — . 
In the phase diagram of CeCoIn 5 depicted in Fig. 1, such 
a crossover is indicated by the hatched area between the 
domains of NFL and LFL behavior and also by the line 

Across (-^0 ■ 

We observe that the end point H c o of the curve 
T CIOSS (H) nominally separating NFL and LFL phases is 
a magnetic-field-induced quantum critical point (QCP) 
hidden in the SC stated. This the most salient feature 
of the phase diagram for the behavior of the resistivity 
p(T,H). Since the entropies of the two phases are dif- 
ferent, the SC transition must become first-order— near 
the QCP, in agreement with the experimental. Moreover, 
LFL behavior remains in effect in the domain T — > 0, 
H > H C 2- It follows that imposition of fields H > H C 2 
will drive the system from the SC phase to the LFL 
phase, where FC, or equivalently the fiat portion of the 
spectrum e(p), is destroyed. Thus, application of a mag- 
netic field H > H C 2 to CeCoIn 5 is predicted to cause a 
step-like drop in its residual resistivity po, as is in fact 
seen experimentally^. Furthermore, it is to expected that 
the higher the quality of the CeCoIns single crystal, the 
greater is the suppression of po- 

As suggested by this analysis of the H — T phase dia- 
gram, the behavior of the dimensionless thermal expan- 
sion coefficient, treated as a function of the dimension- 
less temperature TV, is found to be almost universal. 
Both the left and right panels of Fig. [2] show that all 
the normalized data extracted from measurements on 
CcCoIn j 32 ! 33 collapse onto a single scaling curve. As 
seen from the left panel, the dimensionless coefficient 
C(n{T,H) = a(T 7 H)/a(Tjy 7 H), treated as a function of 
TV = T/T cross , at TV < 1 shows a linear dependence, 
as depicted by the dash-dot line, implying CeCoIns ex- 
hibits LFL behavior in this regime. At TV — 1 the sys- 
tem enters the narrow crossover region. At TV > 1, NFL 
behavior prevails, and both a and S* cease to depend 
significantly on T, with ajq remaining close to the hor- 
izontal line. The observed limiting behaviors, namely 
LFL with a oc T and NFL with a = const., are consis- 
tent with recent experimental results^ 3 ., as it is seen from 
the right panel. From this evidence we conclude that es- 
sential features of the experimental T — H phase diagram 
of CeCoIngi are well represented by Fig. [1] 

In calculations of low-temperature transport properties 
of the normal state of CeCoIns, we employ a two-band 
model, one of which is supposed to be flat, with the dis- 
persion given by Eq. ([3]), while the second band is as- 
sumed to possess a LFL single-particle spectrum having 
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where n is the number density of electrons. 






FIG. 3: Scattering diagrams that contribute to the imaginary 
part of the mass operator E(e), referring to the band with 
finite value of the Fermi velocity. The single line corresponds 
to a quasiparticle of that band, and the double line to a FC 
quasiparticle. 



FIG. 2: (color online). Left panel: Normalized low- 
temperature thermal expansion coefficient otN vs normalized 
temperature Tjv of the normal state of CeCoIns at differ- 
ent magnetic fields H shown in the legend. All the data 
represented by the geometrical symbols are extracted from 
measurements^. The dash-dot line indicates the LFL behav- 
ior taking place at low temperatures under the application of 
magnetic fields. The NFL behavior at higher temperatures 
characterized by both a = const, and S-* of Eq. ((2| is shown 
as the horizontal line. Right panel: Normalized low tempera- 
ture thermal expansion coefficient (a/T) at vs Tn at different 
H shown in the legend. The data^ and T were normalized by 
the values of a/T and by the temperature T„f ~ T crosa , cor- 
respondingly, at the inflection point shown by the arrow. The 
horizontal solid line depicts the LFL behavior, (a/T) = const. 
The dash line displays the NFL behavior {a/T) <x 1/T. 



In conventional clean metals obeying LFL theory, the 
damping j{T) is proportional to T 2 , leading to Eq. ([T]) 
with n = 2. NFL behavior of a{T) is due to the NFL tem- 
perature dependence of j{T) associated with the presence 
of FCSii. In the standard situation where the volume rj 
occupied by FC is rather small, overwhelming contribu- 
tions to the transport come from inelastic scattering, rep- 
resented diagrammatically in Fig. [3^ and Fig. where 
FC quasiparticles (distinguished by the double line) are 
changed into normal quasiparticles; or vice versa, normal 
quasiparticles turn into FC quasiparticles. Contributions 
of these processes to the damping 7 are estimated on the 
basis of a simplified formula^! 



7(p,e)oc J J J J |r(p,Pi,q| 2 ImG H (p - q, e 







finite T- independent dispersion. We begin the analysis 
with the case H = 0, where the resistivity of CeCoIns is 
a linear function of T. As will be seen, this behavior is 
inherent in electron systems having fiat bands. We first 
express the conductivity cr{T) in terms of the imaginary 
part of the polarization operator n(j)2i, 



a = limcj ImII(j, uj — > 0) oc 



1 



dvds 



T J J cosh 2 (e/2T) 
x |r(j,w = 0)| 2 ImG i? (p,e)ImG fl (p,e), (4) 

where dv is an element of momentum space, 7~(j, oS) is the 
vertex part, j is the electric current, and Gr{p, e) is the 
retarded quasiparticle Green function, whose imaginary 
part is given by 



ImG fl (p,e) = - 



1 



e(p)) 2 + 7 2 



(5) 



in terms of the spectrum e(p) and damping 7 referring 
to the band with the finite value vp of the Fermi veloc- 
ity. Invoking gauge invariance, we have T{j,u) = 0) = 
e<9e(p)/9p . Upon inserting this equation into Eq. ((4]) 
and performing some algebra we arrive at the standard 
result 



a{T) 



e 2 n 



7C0' 



(6) 



x ImG_R(— pi, -e)ImGj{(q — pi,u> — e)dpidqdujde, (7) 

where now the volume clement in momentum space in- 
cludes summation over different bands. Calculations 
whose details can be found in£ yield 

7(e) = 77(70 +7ie), ReS(e) = -?77ieln|^j, (8) 

with r\ denoting the volume in momentum space occupied 
by the flat band, and e c being a characteristic constant, 
specifying the logarithmic term in S. Accounting for ver- 
tex corrections^ ensures transparent changes in Eq. ([7| 
and cannot be responsible for the effects discussed in our 
article. We note that Eq. ((5J) leads to the lifetime r q of 
quasiparticles, h/r q ~ ai + a,2T . where h is Planck's con- 
stant, eti and a-i are parameters. This result is in excel- 
lent agreement with experimental observations^. Given 
this result, one finds that p(T) = po + AT, i.e., the re- 
sistivity p{T, H = 0) of systems hosting FC is indeed a 
linear function of T, in agreement with experimental data 
on CeCoIns. Furthermore, the term po arises even if the 
metal has a perfect lattice and no impurities at all. 

The presence of the flat band manifests itself not only 
in kinetics but also in the thermodynamics of CeCoIns, 
e.g., in the occurrence of an additional term AG = C s — 
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C n in the specific heat C(T), given by 



AC = - 



1 

2T 



dA 2 (p) 
dT 



n(p)(l - n(p))du, (9) 



where A (p) is the energy gap and C s and C„ are the spe- 
cific heats of superconducting and normal states, respec- 
tively It is FC contribution that endows CeCoIn 5 with a 
record value of the jump AC /C n = 4.5— of the specific 
heat at T c (the LFL value being 1.43). The enhance- 
ment factor is evaluated by setting T = T c in Eq. ([9]). 
Importantly in systems with flat bands, the quantity 
q = -(l/2T c )(dA 2 /dT) has the same orderi^ as in LFL 
theory, where q ~ 5. To illustrate the point, suppose 
that the momentum distribution n{p) depends only on 
the absolute value of p. One then obtains 



AC(T C ) 
C n (T c ) 



f /-*>(! 



lip)) dp. 



(10) 



Thus, the ratio AC(T c )/C n {T c ), shown to be pro- 
portional to the volume in momentum space occu- 
pied by the flat band, behaves as 1/T C , implying that 
AC(T c )/C n (T c ) diverges at T c — > 0, in agreement with 
data on CeCoIn 5 , where T c is only 2.3 K. 



The formula (GO) can be recast in the forrn^ 



AC(T C ) 
C n (T c ) 



= qS- 1 (T c ) 



TcX(Tc) 
C n (T c ) ' 



(11) 



suitable for extracting the Stoner factor S(T C ) = 
X{T C ) / xo{T c ) ■ Based on the experimental data^, one 
finds S(T C ) ~ 0.3. 

It is straightforward to apply these results to analysis 
of the slope of the peak of the specific heat in the super- 
conducting state of CeCoIn 5 as T — > T c , based primarily 
on Eq. ©. More precisely, one has 



dC(T->T c )/dT 



AC 



dC BCS (T ^T c )/dT C n (T c y 



(12) 



the right side of this relation being evaluated with the aid 
of Eq. (flU)) . This narrowing of the shape of C(T) toward 
the A-point curve, so familiar in the case of superfluid 
4 He, is in agreement with experimental data on CeCoIn5. 

The component of the damping 7 linear in energy 
given by Eq. (|5J) is responsible for a logarithmic cor- 
rection to the specific heat C(T) of the normal state of 
CeCoIn 5 , observed inS&. With f(e) = [1 + e^ 7 ]" 1 and 
R(p,e) = — i In (Gr(p, e)/Ga(p, e)), the formula of^i for 
the entropy S(T — > 0) is recast as 



S(T) = —T~ 



dv 



e^lR(p,e)de. (13) 



where, by virtue of Eq. ©, R(p,e) = 
tan -1 (77/ (1 + rj\n(e c /\e\) — e(p)/e)). Changing vari- 
ables to w = e(p)/e cx (p — Pf)/s and e = zT and 



retaining only leading terms, S(T) is expressed as the 
sum S = S+ + S- , with 



S+ oc T 



(l + e z ) 2 



tan ; — -. — 7—7 dw. 



l + V \n(e c /T)Tw 



(14) 

These integrals are evaluated analytically to yield S(T) — 
Sfl{T) oc rjThiT, in agreement with available experi- 
mental data on the specific heat of CeCoIng^. 

Now we show that the application of field H > H C 2 
on CeCoIns generates the step-like drop in the residual 
resistivity pq. Indeed, as seen from Fig. [U at low tem- 
peratures T < T cross , the application of fields H > H C 2 
drives the system from the SC state to the LFL one, 
where the flat portion of e(p) is destroyed. Thus, the 
term rj vanishes strongly reducing pq. The suppression 
of po in magnetic fields is associated with another NFL 
phenomenon recently observed in CeCoIn^. This is the 
deviation of the temperature dependence of the ampli- 
tude of magnetic oscillations A(T) from the standard 
Lifshitz-Kosevich-Dingle form 



A(T) = A D (X/sinhX) 



(15) 



Here Ad oc e Y is the Dingle factor, while Y = 
277 2 k B T D /u} c , X = 2n 2 k B T/uj c , T D is the Dingle tem- 
perature, and u) is the cyclotron frequency Convention- 
ally Tq, assumed to be constant, is associated with im- 
purity scattering. More generally, Tjj is related to the 
value of po in Eq. ((T|), which, as we have seen, becomes 
T- dependent in a system with a flat band. Why? Be- 
cause the flat band, in the present context, guarantees 
an overwhelming contribution to po that must disappear 
at T < T CYOSS (H) to evade violation of the Nernst theo- 
rem. This inexorable sequence triggers the abrupt down- 
ward jump of po and of T B correspondingly, leading in 
turn to an upward jump of the Dingle factor Ad near 
Tcross (H), which agrees with observations^. In parallel 
with this challenging behavior of the residual resistivity 
in magnetic fields, it is apposite to address the behavior 
of po versus pressure P, studied in CeCoIns experimen- 
tally in^. At P > P* =1.6 GPa, po drops reversibly by 
one order of magnitude to a very small value of about 
0.2 /iOcm. We can reasonably infer that the application 
of pressure eliminates FCS, triggering a jump in po to 
the lower value measured. It should be emphasized that 
a nonzero contribution of FC to po is associated with the 
presence of the crystal lattice, more precisely, with the 
Umklapp processes, violating momentum conservation. 
At the same time, such a restriction is absent in dealing 
with the thermal resistivity wq- If, as usual, one normal- 
izes the thermal resistivity by w = ir 2 T/(3e 2 K,) where k is 
the thermal conductivity, the famous Wiedemann- Franz 
relation then reads po = Wq. The distinguished role of 
the Umklapp processes in the occurrence of po in Fermi 
systems with FC implies that in the presence of FC, the 
Wiedemann-Franz law is violated so that po < Wq. In 
CeCoIns this violation does take the placed. 
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In summary, we have shown that the application of 
magnetic fields and pressure on CeCoIn 5 leads to strong 
suppression of the residual resistivity p . By considering 
the behavior of the thermal expansion coefficient, the spe- 
cific heat, and the amplitude of magnetic oscillations, we 
have unveiled the roles played by the flat band in thermo- 
dynamic as well as in transport properties of CeCoIn5. 
Our considerations furnish strong evidence for the pres- 



ence of flat band in CeCoIns , which thereby becomes the 
member of a long-expected class of Fermi liquids. 
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